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Abstract 

A systematic correlation between the initial profile of discrete breathers and their 
frequency is described. The context is that of a very weakly harmonically coupled chain 
of softly anharmonic oscillators. The results are structurally stable, that is, robust under 
changes of the on-site potential and are illustrated numerically for several standard 
choices. A precise genericity theorem for the results is proved. 

1 Introduction 

Chains of coupled anharmonic oscillators have many applications in condensed matter and 
biophysics as simple one-dimensional models of crystals or biomolecules. Of particular 
interest are the so-called "breather" solutions supported by such chains, that is, oscillatory 
solutions which are periodic in time and exponentially localized in space [0. The simplest 
possible class of models, where all the oscillators are identical, with one degree of freedom, 
and nearest neighbours are coupled by identical Hooke's law springs has been widely studied 
in this context. The equation of motion for the position q n (t) of the n-th oscillator (n G Z) 
is 

q n ~ a(q n+1 - 2q n + q n _{) + V'{q n ) = (1) 

where a is the spring constant and V is the anharmonic substrate potential, which we 
choose to normalize so that V'(0) = and V"(0) = 1. MacKay and Aubry have proved 
the existence of breathers in the weak coupling (small a) regime of this system [^]. They 
noted that in the limit a — > 0, system (|l]) supports a one-site breather, call it q = (q n )nez, 
where one site (n = say) oscillates with period T while all the others remain stationary at 
0. Using an implicit function theorem argument they proved existence of a constant period 
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continuation q a of this breather away from a = provided a remains sufficiently small and 
T £ 2vrZ+. 

Having established the existence of these breathers, the question remains: what do 
they look like? Of course, this problem has been intensively studied, and several efficient 
numerical schemes have been devised [||, ||]. The purpose of this note is to point out that 
significant, systematic qualitative information about breathers in weakly coupled chains can 
be obtained from the initial direction of continuation of one-site breathers, 

, <9q Q 

qo - 



9a a=0 

Specifically, the small a behaviour of breathers can be determined by approximating the 
curve q a by its tangent line at a = 0, that is 

q« = qo + "Qo + o(a). (3) 

Note that equation (|3|) is a precise mathematical statement, following directly from the 
MacKay-Aubry theorem, not a formal expansion. This allows one to make a systematic 
study of the frequency dependence of breather initial profiles for several on-site potentials. 
The results show a generic behaviour of alternating frequency bands of two qualitatively 
different types of breather, which we call "in-phase breathers" and "anti-phase breathers" . 
We go on to prove that this alternating behaviour is generic in a precise sense. 



2 The direction of continuation of one-site breathers 

In the following we will assume that V : R — ► R is twice continuously differentiable and has 
a normalized stable equilibrium point at (V(0) = 0, V"(0) = 1). Consider the equation 
of motion for a particle moving in such a potential, 

x + V'{x) = (4) 

with x(0) = 0. Provided |x(0)| is small enough, x(t) must be a periodic oscillation. All the 
potentials we consider will be softly anharmonic with classical frequency spectrum (0, 1). 

Let xx(t) denote the solution of @ with period T > 2tt and x(0) > which has even 
time-reversal symmetry, xt(— t) = XT(t). From this we may construct a 1-site breather 
solution of system (||) with a = 0, call it qo: 

U s ( x T (t) n = . . 

M*) = ( o n/0 . (5) 

The MacKay-Aubry theorem establishes the existence of a continuation q Q of this solution 
away from a = in a suitable function space, defined as follows. 

Definition 1 For any n E N, let C^ n denote the space of n times continuously differen- 
tiable mappings R — ► M which are T -periodic and have even time reversal symmetry. Note 
that C^ n is a Banach space when equipped with the uniform C n norm: 

\q\ n :=sup{|g(i)|,|g(i)|,...,k (n) (i)|}. 
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Definition 2 For any n G N, 

:= {q : Z — > Cy n such that ||q|| n < 00} 

u>/iere 

||q|| n := sup \q m \ n . 

mEZ 

iVote i/iai (^<£ n , || • \ \ n ) is also a Banach space. 

The required function space is J7^ 2 . By construction, qo £ n^ 2 and is exponentially 
spatially localized. 

Theorem 3 (MacKay-Aubry) If T £ 27rZ f/iere exists e > swc/i i/ta£ /or a// a £ [0, e) 
i/iere is a unique continuous family q a G fij 2 °/ solutions of system |Jj at coupling a 
with qo as defined in ffl). These solutions are exponentially localized in space and the map 
[0, e) — > Q^ 2 c/uien 6y a i— > q Q is C 1 . 

The idea of the proof is to define a C 1 mapping : 0^ 2 © R — > 0^ , 

F(q, a) m = g" m - a(g m +i - 2g m + g m -i) + V(g m ), (6) 

so that .F(q, a) = if and only if q is an even T-periodic solution of system (|l|). In 
particular, ^(qo, 0) = by construction. Using T (jz 2-7rZ + and anharmonicity of V, one can 
show that the partial derivative of F with respect to q at (qo,0), DF qo : £1^2 ~ * ^to> * s 
invertible (DF qo is injective with (DF qo )~ l bounded). Hence the implicit function theorem 
H applies and local existence and uniqueness of the C 1 family q Q satisfying 

F(q a ,a)=0 (7) 

are assured. Persistence of exponential localization is proved as a separate step. 

The object of interest in this paper is q G f2y 2 > the tangent vector to the curve q a 
at a = 0. This may be constructed by implicit differentiation of (|?]) with respect to a at 
a = 0: 

_! dF 
da 



dF 



q = 



(8) 

(qo,o) 



(qo,0) 

Both -Di^o and dF/dalr^ q\ are easily computed from (|5]||). Evaluating the right hand side 
of (||) to find q = % G * s then equivalent to solving the following infinite decoupled 
set of ODEs for {xm £ Cy 2 : m G Z}: 

Xo + ^"(xt)xo = -2x T (9) 

Xm + Xm = X T \m\ = 1 (10) 

Xm + Xm = \m\ > 1. (11) 

Recalling that T ^ 27rZ + (so cos ^ C^t 2 ) one sees horn ( |TT] ) that x m = for |m| > 1. So one 
need only solve the linear boundary value problems (PJl0|) with x m (0) = 0, Xm(T) = Xm(0) 
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for m = 0, 1 (clearly x-i = Xi)i coupled to the nonlinear driver equation (|j) for xt- This 
is a numerically trivial task. 

Having computed the initial value of the tangent vector, we may approximate the 
breather initial profile for small a using (||): 

{x T (0) + axo(0) + o(a) m = 
«Xi(0) + o(a) |m| = 1 (12) 
o(a) |m| > 1. 

So to first order in a, the continuation leaves all but the central (m = 0) and off-central 
(m = ±1) sites at the equilibrium position. The qualitative shape of the breather initial 
profile depends crucially on xi(0) (but not Xo(0))- If Xi(0) > 0, the continuation displaces 
the off-central sites from equilibrium in the same direction as the central site. The result is a 
hump shaped breather in which the central and off-central sites oscillate, roughly speaking, 
in phase (that is they attain their maxima and minima simultaneously). We shall call such 
breathers "in-phase breathers (IPBs)." If Xi(0) < 0, on the other hand, the off-central sites 
are displaced in the opposite direction from the central site, resulting in a sombrero shaped 
initial profile. In this case, the central and off-central sites oscillate, roughly speaking, 
in anti-phase (the central site attains its maximum when the off-central sites attain their 
minima, and vice- versa). We shall call such breathers "anti-phase breathers (APBs)." We 
shall find that for a generic substrate potential V, system (|l|) supports both IPBs and APBs, 
the type varying in bands as T increases through (2ir, oo). 



3 Numerical results 

The results presented in this section were generated using a 4th order Runge-Kutta method 
with fixed time step 5t = 0.01 to solve the driver equation (|4|) for xt and the linear boundary 
value problems (PJTo|). The results obtained depend crucially on whether the substrate 
potential has reflexion symmetry about the equilibrium position. 

The following phenomenologically standard asymmetric potentials were investigated: 

Morse: V M {x) = -(1 - e~ x ) 2 (13) 



1 / 1 2 

Lennard-Jones: V LJ {x) = ^[^n ~ ~&) ( 14 ) 



Cubic: V c {x) = ^x 2 - ^x 3 . (15) 

Note that, although even, Vlj is asymmetric about its equilibrium position x = 1. Graphs 
of Xo(0) and xi(0) against T for all these potentials are presented in figure |l[ 

Since xo is always positive, continuation always proceeds initially by pulling the central 
oscillator further away from equilibrium as a increases. Note also that xo remains bounded 
as T — > oo if the single oscillator solutions xx{t) remain bounded (Vc) but grows unbounded 
otherwise (Vm and Vlj). In all cases it grows unbounded as T — > 2tt. 

From figure [l] we see that the sign of Xi(0) changes as T increases, leading to the 
prediction of T-bands of IPBs and APBs, as explained in section 0. The sign changes 
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Figure 1: Graphs of Xo(0) (left) and xi(0) (right) against T for various asymmetric on-site 
potentials (solid: Morse; dashed: Lennard- Jones; dotted: cubic). 



are clearly associated with the vertical asymptotes of the graphs at each T S 27rZ + . The 
presence of such asymptotes may be understood by using Green's function techniques to 
write down First note that the particular integral for ( jlOj ) with initial data z(0) = 

i(0) = is 



z{t) 



I sin(i — s)xt(s) ds. 
Jo 



(16) 



It follows from evenness of \i that Xl(*) = Xi(0) cos * + z (t)i anci hence, applying the 
periodicity condition xi(0) = Xi(T), that 



Xi(0) 



1 f T 

/ cos(T — t)xT{t) dt. 

sin T ,/ 



(17) 



Equation (|17|) shows that there is a vertical asymptote at T = 2mr, with a sign change in 
Xi(0), unless 



r2nir 

/ COSt X2nir(t) dt = 0, 

Jo 



(18) 



that is, unless X2n-rr{t) has vanishing n-th Fourier coefficient. It might appear from (|T^) 
that there should also be asymptotes at T = (2n + l)ir. This cannot be true since standard 
results on continuity of solutions of ODEs with respect to initial data imply that xi (0) is 
continuous for T ^ 2-7rZ + . In fact, a simple argument using periodicity and evenness of xt 
demonstrates that 

-(2n+l)7T 

COS t X( 2n +l)n(t) dt = (19) 



for all n G Z + and V . By contrast, we will prove in section |^ that (|18|) almost never holds 
(in a sense which will be made precise), so that the resonant periods T = 2mr generically 
separate IPB bands from APB bands. 

Turning to symmetric potentials, the following standard examples were investigated: 



Frenkel-Kontorova: Vfk{ x ) =1 — cosx 



(20) 



Figure 2: Graphs of Xi(0) against T for various symmetric on-site potentials (solid: Frenkel- 
Kontorova; dashed: quartic; dotted: Gaussian). 



Quartic: Vq{x) 
Gaussian: Vg(x) 



1 2 
-X 

2 

1 - 



1 4 

-r 

-x 2 /2 



(21) 
(22) 



The graphs of xo display identical behaviour to that found for the asymmetric examples, and 
hence are not illustrated. The graphs of Xi(0) against T, on the other hand, are strikingly 
different from the asymmetric case (see figure |^). In each case, vertical asymptotes are 
present at T = 2mr only if n is an even positive integer. To explain this, note that evenness 
of V implies that XT(t) is T/2 antiperiodic, that is xx(t — T/2) = —XT(t), which in turn 
guarantees that equation (|l8|) holds whenever n is odd. So alternating bands of IPBs and 
APBs occur for these potentials, but the alternation occurs half as often. 

One universal feature of the graphs of xi (independent of symmetry) is that xi ~* 00 
as T — > 2tt. This may easily be understood in the case where V is analytic (as in all our 
examples) by constructing a Linstedt expansion for xt(0 in the small amplitude regime Q. 



4 Genericity of numerical results 

The aim of this section is to prove the generic existence of a sign change in xi as T crosses 
a resonance T G 2ir7* + . In particular, we will prove that, for a given n, the set of potentials 
on which (|l8|) holds is negligibly small. To be precise, we will show that, for potentials in a 



certain Banach space V, the subset of potentials on which (18) holds is locally a codimension 
1 submanifold. 



Definition 4 Let V = {V : R - 

be the mapping 

F(q,V)=q + V'(q). 
Note that F is a C 1 mapping between Banach spaces. 



such that V is C 2 and \\V\\2 < oo} and F : C^ 2 ®F 
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Definition 5 A potential V £ V is anharmonic at q £ Cj^ 2 if Q * s nonconstant, F(q,V) = 
0, and kei DF( q y) = {0}, where DFi q y\ : C^ 2 ~^ Cro denotes the partial derivative of F 
at{q,V). 

Clearly F(q, V) = means that q is a solution of Newton's equation for motion in poten- 
tial V. That injectivity of DFr q y\ is equivalent to the standard definition of anharmonicity 
in Hamiltonian mechanics is shown, for example, in the original proof of Theorem |3] [Q]. For 
all the potentials we considered, V is anharmonic at q = xt £ C*t 2 ^ or a ^ ^ ^ 7r - 

Definition 6 ^4 perturbation neighbourhood of (q, V) £ „ (B V is a pair (hi, f) where 

U C V is an open set containing V , and f : IA — > C^ 2 *s a C 1 map satisfying (a) every 
W £ U is anharmonic at f(W) and (b) fiV) = q. Note that U is trivially a Banach 
manifold. 

Lemma 7 For any V £ V anharmonic atq<E C^ 2 there exists a perturbation neighbourhood 
(U, f) of (q, V). The function f is unique on sufficiently small U. 

Proof: Anharmonicity implies that DF^ q y^ is injective, and hence the standard solvability 
criterion in linear ODE theory guarantees it is also surjective (see section 3.3.2 of || for 
details). The open mapping theorem then ensures boundedness of DF^yy Hence DF( q V ^ 
is invertible, and the implicit function theorem applied to F ensures existence of (Li, f) and 
local uniqueness of /. □ 

Definition 8 For T £ 2-7rZ + ; we shall say that V £ V is degenerate at q £ 2 if V is 
anharmonic at q and 

T 

cos t q(t) dt = 0. 

Theorem 9 Let T £ 27rZ + , Vq £ V be degenerate at qo £ 2 , and {Li, f) be a perturbation 
neighbourhood of (qo,Vo). The subset oflA on which degeneracy persists is a codimension 1 
submanifold. 

Proof: Consider the C 1 mapping / : U — > R defined by 

I(V) = f T cost [f(V)] (t)dt 
Jo 

(I is C 1 since / is C 1 ). By the Regular Value Theorem (see appendix) the result follows 
if we establish that is a regular value of /, or in other words that I is a submersion at 
every V £ 7 -1 (0). Since ker DLy is of finite codimension in TylA = V it splits and hence it 
suffices to show that Dly : V — > K is surjective for all V £ / -1 (0). 

Let V £ I' 1 ^), Q = f(V). For any SV £ V, let Sq sv = -[DF^ V) (SV' o q)} £ C+ 2 . In 
other words 5qsy is the unique even, T-periodic C 2 solution of 

ksv + V"(q(t))8q 5V = -5V'(q(t)), 
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whose existence follows from anharmonicity of V at q (invertibility of DFr^ y^ on C^ ). 
With this notation, 

DI V (6V)= / costSq 6V (t)dt 
Jo 

and surjectivity of Dly will follow if we exhibit 5V G V such that DIy(5V) ^ 0. We 
construct such a <5V in two steps. 

For some to G (0, T/2) and e > 0, let 6 G C^ 2 be a non-negative function with supp&n 
[0, T/2] = [t - e, t + e]. Clearly 

f T f T/2 

/ 6(t) cos t dt = 2 / 6(t) cos tdt^O 
Jo Jo 

(where T G 2-7rZ + has been used) provided we choose to ^ cos _1 (0) and e sufficiently small. 
Since the critical points of q{t) are isolated, we may also assume that q is invertible on 
[to — e, to + e] with C 2 inverse q -1 : \x\ , 22] — > [to — e, to + e]. By its construction the function 
<7 defined by 

„M _ / -[P F (9,y) b ] °q~ X ]{x) x G [si,x 2 ] 
51 ' 1 x^[xi,x 2 ] 

is in C°(R) (since DF maps to Cj, ). This allows us to construct a function 8V{x) G C 1 (IR) 
such that Sq$ v = b as follows. For x G ran(q), 5V"(x) = g(z) dz defines a C 1 function on 
ran(g). Since Sq sv is clearly independent of the behaviour of 5V outside ran(g), we extend 
SV to a C 1 function on R with support contained in some compact set K D ran(g). 

We cannot immediately conclude that Dly is surjective, since 5V is C 1 but not nec- 
essarily C 2 . However, using a density argument it is straightforward to find a SV G C 2 
close enough to SV such that DIy(5V) 7^ stills holds. More precisely, Dly has a natural 
continuous extension, call it Dly to C . By construction 6V satisfies DIy(5V) 7^ and 
belongs to C\ (the C 1 functions with support contained in K). Since C|- is dense in Cj^, 
there exists a sequence 5V m G C\ C "P converging (in C 1 norm) to and by continuity 
DIy(5V m ) 7^ for all m sufficiently large. □ 

In the above, we have chosen V = (C 2 (R),|| • [[2) as our space of potentials. This 
choice was made for the sake of clarity and notational simplicity - many other choices 
would work. In fact, all but two (Vfk an d Vg) of the example potentials considered in 
section || lie, strictly speaking, outside V, since they are unbounded. However, the analysis 
above can easily be adapted to deal with this: one simply replaces V by the affine space 
Vy := {W : \\W - V\\ 2 < 00}. 



5 Comparison with tail behaviour 

All our considerations so far have concerned the frequency dependence of breather initial 
profiles in the breather core, \m\ < 1. It is interesting to compare these results with the 
profile behaviour in their tails, that is, for large \m\. Theorem |3] guarantees that breathers 
decay spatially exponentially fast, but gives no explicit information on the decay rate, or 
on sign correlations between neighbouring sites in the tail. Such information is obtainable, 
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at least heuristically, by means of the Fourier analytic approach of Flach M. We briefly 
recall this computation, adapted to our specific context. 

Since the breather is spatially localized, time periodic with period T = 2ir/u and time 
reflexion symmetric, it has a Fourier series of the form 



oo 



(t) = y A m ^ cos hut, (23) 



k=0 



with A m k — > as m —* ±00 for all k. Substituting (||) into ([[]) one obtains a nonlinear 
algebraic system for the coefficients A m ^. Given the spatial decay criteria one expects these 
to be close, for large \m\ to solutions of the linearized algebraic system. This linearization 
decouples into a separate 2nd order linear difference equation for each mode k, namely 



A 



m+l,k 



2 + 



l - kW 



A m , k + A m _ hk = 0, (24) 

a 

whose general solution lies in the span of /^T m where \x k is a root of 

1 - k 2 u? 

A* 2 - 6^ + 1 = 0, £k ■= 2 + (25) 

a 

with \fik\ < 1- The condition A m ^ — > implies [ik is real and \/j,k\ < 1- Hence the coefficient 
£fc in ( |25| ) must lie outside [—2, 2] for all k, so kuj ^ [1, \/l + 4a]: every harmonic of to must 
lie outside the phonon frequency band. Given this, all Fourier modes decay exponentially 
like ^ ' m '- 

In order to understand the tail behaviour, one must identify the dominant harmonic 
for large \m\, that is, the harmonic with slowest decay. Algebraically we seek that k, call 
it k^ say, which has |/Ufc| closest to 1. Thinking of ^ as the preimage / -1 (£fc) where 
/ : (—1, 1)\{0} — > R, f(fj,) = fi + fj, , one sees that this is the k for which |^| is closest to 
2. If > 2 (kuj < 1, below the phonon band), then 

a|& - 2| = 1 - (kuo) 2 , (26) 



while if £fc < — 2 {kuj > vT + 4a 2 , above the phonon band), then 

a|e fc + 2| = (l + 4a)-(M 2 - (27) 

Hence, the dominant harmonic is the one lying closest to the phonon band in squared 
frequency space. If kdOJ lies below the band, then < /ife d < 1, and the tail has spatially 
uniform sign - the analogue for tail behaviour of an IPB. If k^co lies above the band, then 
— 1 < fik d < so the breather has an alternating tail, the tail analogue of an APB. 

It is straightforward to partition the breather existence domain in (T, a) parameter 
space into IPB and APB subdomains according to this analysis, as shown in figure ^. In 
the period window T £ [2nk, 2ir(k + 1)], one site breathers may be continued vertically from 
a = at most until the phonon band captures the {k + l)th harmonic, (fc + \)oj = \/l + 4a, 
giving the upper bounding curve 

T ex (a) = 2tt^^=. (28) 
yl + 4a 
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T/n 

Figure 3: The breather existence domain for generic choice of substrate, partitioned into 
IPB-tail and APB-tail sectors, denoted + and — respectively. 

This is the jagged upper curve depicted in figure |J| Within the existence domain below this 
curve, the transition between IPB and APB occurs where the pair kuj, (k + l)u> straddling 
the phonon band become equidistant (in squared frequency space) from that band: 

\-k 2 u? = (k + 1) 2 lu 2 - (1 + 4a) 

^T trans {a) = 2vry - — . (29) 

These boundaries are depicted as dashed curves in figure ||. For fixed a, the transition 
with increasing T across each boundary is always from IPB to APB, denoted + and — 
respectively in the diagram. 

Comparing this picture with figure |l| (which should be thought of as valid in a thin 
strip of the (T, a) plane containing the T axis) one sees that there will be many period 
intervals in which the core IPB/APB classification conflicts with the tail classification. In 
defence of our choice of terminology, based on the core behaviour, we note that breathers are 
interesting precisely because they are localized modes whose tails are exponentially small, 
so the core behaviour should be regarded as more significant than the tail behaviour. It 
is also worth noting two similarities between the core and tail analyses at small a. First, 
they are structurally robust, that is, qualitatively invariant under generic changes of the 
substrate potential. Second, the resonant frequencies to = \/n mark transitions between 
IPB and APB behaviour in both the core and tail senses. 

Several caveats should be attached to the picture developed in this section. First, 
theorem || gives no guarantee that breathers really will exist in the whole domain depicted. 
Generically, it appears that they do, but there are Klein-Gordon-like systems with smaller 
existence domains ||. Second, we have predicted the decay properties of a breather's Fourier 
modes based on the linearized model. Numerically this approximation seems to work well 
most of the time. However, occasionally a few modes have anomalously slow spatial decay 
due to nonlinear corrections ||], so figure || may receive fine detail corrections. Third, it 
may be that some of the Fourier modes are are wholly absent for symmetry reasons (e.g. 
V is symmetric). The picture is then modified substantially, in much the same way as the 
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core picture changes from figure [T] to figure 0. Finally, these results, while plausible, are 
heuristic. It would be interesting to seek a rigorous analytic underpinning for them. 

6 Concluding remarks 

In this paper, a systematic correlation between breather initial profiles and their frequencies 
was examined. Two types of breather were identified, called IPBs and APBs, based on sign 
correlations of the 3 core lattice sites. The numerical data suggest that generically these 
two types occur in alternating bands in the T parameter space (T £ (2ir, oo)), and that the 
resonant periods T S 27rZ + separate an IPB band from an APB band. The genericity of 
this behaviour was proved rigorously. 

The distinction between IPBs and APBs, which does not appear to have attracted 
much attention in the literature, may have phenomenological implications in applications 
of the model (Q). One reason for interest in discrete breathers (particularly continued one- 
site breathers) is that, because of their strong spatial localization, they typically require 
little energy to achieve large amplitude oscillations close to the centre. This makes them 
good candidates for "seeds" of mechanical breakdown of the network, the idea being that 
the central oscillation becomes so violent as to break the chain (a similar mechanism is 
postulated as a mechanism for DNA denaturation, for example ||). One expects the central 
intersite springs to carry a larger proportion of the total breather energy for APBs than 
IPBs, making them more energy efficient seeds of mechanical breakdown. However, further 
detailed simulations of full lattice systems are required to test this hypothesis. 
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Appendix 

We recall some basic definitions and an elementary result of infinite dimensional differential 
topology, the Regular Value Theorem, which gives conditions under which the level set of 
a smooth function is guaranteed to be a manifold. 

Definition 10 We say that a closed subspace S of a complete topological vector space B 
splits if there exists another closed subspace C which is complementary to S, i.e. C + S = B 
and CHS = (0). 

Any finite dimensional subspace (necessarily closed) of a Banach space splits, as does any 
closed subspace of finite codimension. 
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Definition 11 AC 1 map f : X — > Y between Banach manifolds is a submersion at x G X 
if Df x : T X X — > Tf( x ^Y is surjective and ker Df x splits. A value y £ Y is a regular value 
of f if f is a submersion for every x £ f^ 1 (y)- 

Note that in the case that Y is finite dimensional then ker Df x is of finite codimension in 
T X X and hence is guaranteed to split. In this case we need only verify that Df x is surjective 
for / to be a submersion at x. 

Theorem 12 (The Regular Value Theorem) Let f : X — > Y be a C 1 map between Banach 
manifolds and y G Y be a regular value of f. Then is a submanifold of X with 

T x f- 1 (y)^ker Df x . 

The proof follows from working in charts around x and f{x) and applying the Implicit 
Function Theorem |1C]. 



References 

[1] S. Flach and C.R. Willis, "Discrete breathers" Phys. Rep. 295 (1998) 181. 

[2] R.S. MacKay and S. Aubry, "Proof of existence of breathers for time- Reversible or 
Hamiltonian networks of weakly coupled oscillators" Nonlinearity 7 (1994) 1623. 

[3] J.L. Marin and S. Aubry, "Breathers in nonlinear lattices: numerical calculation from 
the anticontinuous limit" Nonlinearity 9 (1996) 1501. 

[4] T. Bountis, H.W. Capel, M. Kollmann, J.C. Ross, J.M. Bergamin, and J. P. van 
der Weele, "Multibreathers and homoclinic orbits in 1-dimensional nonlinear lattices" 
Phys. Lett. 268A (2000) 50. 

[5] M. Haskins and J.M. Speight, "Breathers in the weakly coupled topological discrete 
sine-Gordon system" Nonlinearity 11 (1998) 1651. 

[6] Y. Choquet-Bruhat, C. DeWitt-Morette and M. Dillard-Bleick, Analysis, Manifolds 
and Physics, Part I (North-Holland, Amsterdam, 1982) p 91. 

[7] R. Grimshaw, Nonlinear ordinary differential equations (Blackwell Scientific Publica- 
tions Ltd., Oxford, 1990) Chapter 7. 

[8] S. Flach, "Obtaining breathers in nonlinear Hamiltonian lattices" Phys. Rev. E51 
(1995) 3579. 

[9] E.W. Prohofsky, K.C. Lu, L.L. Van Zandt and B.F. Putnam, "Breathing modes and 
induced resonant melting of the double helix" Phys. Lett. 70A (1979) 492. 

[10] Y. Choquet-Bruhat, C. DeWitt-Morette and M. Dillard-Bleick, op. cit. p550. 



12 



